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We show how the configuration-space form of the Bogoliubov ground state wave function of a 
bosonic condensate with a single vortex in a harmonic trap can be described in terms of bosonic 
Jastrow correlations. We then generalize this result to study the first effects of such correlations 
on a mean-field vortex lattice state and show that the included correlations lower the energy below 
that of the mean-field state. Although the reduction is relatively small, it is a precursor of the more 
general expected effect of correlations in describing the melting of the vortex lattice at large angular 
momentum per particle. 



I. INTRODUCTION 



With increasing rotation, the ground state of an ultra- 
cold gas of bosons in a harmonic trap undergoes a tran- 
sition from a vortex lattice with broken rotational sym- 
metry, for which mean-field theory provides a good de- 
scription (see Ref. [l[ and references therein), to a series 
of symmetry-restored and strongly correlated states 
bosonic analogs of quantum Hall states Q . This transi- 
tion is mediated by the correlations present in the inter- 
acting system which are absent in mean-field theory. To- 
wards understanding how these correlations lead to more 
favorable states, we studied in Ref. [IJ a condensate with 
a single vortex; including Bogoliubov fluctuations around 
the mean-field ground state, we showed that the corre- 
lations induced by these fluctuations lower the energy 
of the Bogoliubov ground state compared to that of the 
mean-field ground state and cause an uncertainty in the 
position of the vortex. 

Here we take a first step in generalizing this earlier 
result to a vortex lattice, showing how correlating two 
particles in the wave function lowers the energy of the 
lattice. We focus on correlations described by simple 
bosonic Jastrow factors, (z,i — Zj) 2 , in the wave function, 
where z ~ x + iy is the position of a particle in the 
complex plane. Such correlations tend to lower the in- 
teraction energy by keeping the particles apart and are, 
hence, favored by repulsive interactions. On the other 
hand, each factor carries two units of angular momen- 
tum and, therefore, tends to increase the kinetic energy 
of the system. With increasing angular momentum, more 
and more of these factors enter the wave function, and 
the states become more strongly correlated, e.g., as in 
the Read-Rezayi [(| and the bosonic Laughlin [f| states. 
The evolution of the system as its angular momentum 
increases towards and beyond the melting transition and 
the role that Jastrow correlations play in this phase tran- 
sition is still an open problem [l| . With increasing angu- 
lar momentum, particles begin to occupy single-particle 
states which previously were empty (or had vanishingly 
small occupations) ; this increase in the size of the config- 
uration space of the particles can lead to possible quasi- 
degeneracies between states with different single-particle 
occupations and, consequently, to the onset of quantum 



fluctuations which ultimately destroy the vortex lattice. 
Even for low angular momenta, where only a few vortices 
are present, correlations still play a significant role in re- 
distributing the particles among single-particle states, as 
Cremon et al. find 7] by studying few-vortex (N v < 4) 
systems numerically and comparing the exact and mean- 
field ground states. This paper sheds a complementary 
light on such redistributions. 

An informative example of how correlations function 
is the gas of attractive bosons studied in Ref. @. For 
a total (arbitrary) angular momentum hL, the ground 
state wave function is 



ipi(z;L) = z l c 



(1) 



where z c = X)i=i z %l^ is the center of mass coordinate. 
This state has a total interaction energy ~ —N(N—l)/2. 
Moreover, the wave function 



i' 2 (z; L) = ^2 ( z n ~ z ii) 2 ^i( z ; L-2) 



(2) 



»1<»2 



describes an excited state of the system with the same 
angular momentum and with a higher interaction energy 
^ — N(N — 2)/2. If we now change the nature of the in- 
teractions from attractive to repulsive, these two states 
switch places in the energy spectrum, with ^(L)) be- 
coming lower in energy than \ipi(L)), although it does 
not become the ground state. Similar to the Bogoliubov 
single- vortex state, the Jastrow correlations included in 
\tp2(L)) help to lower the now repulsive interaction en- 
ergy. 

In this paper, we first show how the real-space form of 
the Bogoliubov ground state of the single- vortex conden- 
sate [4[ includes two-particle bosonic Jastrow factors and 
can be expanded as a sum over symmetric polynomials 
with successive number of Jastrow factors. The first term 
of the sum is just the original uncorrelated mean-field 
wave function, and the last term has N/2 simultaneous 
Jastrow factors. The effect of such Jastrow correlations 
is to reduce the total energy by a term 0(N~ 1 ) which, al- 
though small, is a precursor of the more general expected 
effect of correlations. We then generalize the correlated 
single-vortex case to a vortex lattice, initially described 
as a mean-field condensate. Again, we find that the in- 
cluded Jastrow correlations lead to a relative reduction 
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of the energy 0(N~ 1 ). We also find that the inclusion 
of these correlations in the trial wave function leads to a 
nonvanishing density at the vortex cores, indicating the 
presence of quantum fluctuations of the vortices, simi- 
lar to the case of the single- vortex system we previously 
studied 0]. 

In the next section, we delineate the basic model de- 
scribing a condensate in terms of Landau levels. In 
Sec. lIIIl we expand the Bogoliubov ground state in terms 
of a series of TV-particle Fock states with increasing num- 
ber of particles in the two single-particle states connected 
to the mean-field ground state through the interactions; 
we show that these Fock states are represented by mono- 
mial symmetric polynomials (which, in turn, can be ex- 
panded in terms of other symmetric polynomials contain- 
ing Jastrow factors) and find the form of the correlations 
present in the wave function. In Sec. IIV1 we general- 
ize this construction to a vortex lattice system. Finally, 
in Appendix we derive a general algebraic identity 
connecting the monomial symmetric polynomials encoun- 
tered in this problem to symmetric polynomials with suc- 
cessive number of Jastrow factors, and in Appendix [B] 
we lay out the details of the derivations used to arrive at 
the results of Sec. IIIII 



II. BASIC MODEL 

We consider a gas of N bosons of mass to in a har- 
monic trap of frequencies wj_ in the x—y plane and lu z in 
the z direction, rotating around the z axis with angular 
velocity fi. We assume weak two-body repulsive inter- 
actions of strength g = Aith 2 a/m, where a is the s-wave 
scattering length. The Hamiltonian in the rotating frame 
is thus 



N 

E 

i=l 



21 

2m 



1 



+ -m(w^a; 2 + ui\y 2 + oj 2 z 2 ) - Mi 



i<j 



(3) 



where I = z • (r x p) is the angular momentum along the 
z direction. 

In the limit of fast rotation (f2 < at zero tem- 
perature, the gas becomes quasi-two-dimensional and re- 
sides in the axial ground state of the harmonic trap. 
The single-particle eigenstates of the non-interacting sys- 
tem are the Landau levels, |nm), where n is the radial 
quantum number and to > — n is the angular momen- 
tum along the rotation axis. The characteristic interac- 
tion energy scale is Vq — <?/[(27r) 3 / 2 d^d 2 ] where d±^ z = 
yJh/muji_ tZ are the characteristic oscillator lengths in the 
transverse and axial directions. We assume the interac- 
tions to be sufficiently weak that Vb <C 2Hlu±; there- 
fore, as — ¥ ljx, the system resides in the lowest-energy 
(n = 0) manifold of Landau levels. The wave function 
of a particle in the lowest Landau level (LLL) with m 



units of angular momentum, corresponding to the single- 
particle state |0m), is 



(f> m (z) = (z\0m) = 



1 



-M72 



(4) 



where z = (x+iy)/d± is the dimensionless position in the 
complex plane. For brevity, we suppress, throughout this 
paper, the factor exp [—J2iLi |^ | 2 /2]/(d_L- v /7r) Ar common 
to all A^-particle LLL wave functions. 



III. CORRELATIONS IN THE SINGLE- VORTEX 
BOGOLIUBOV GROUND STATE 

We studied in Ref. [|| the properties of a single- vortex 
system in the LLL by including small-amplitude Bogoli- 
ubov fluctuations about a mean-field condensate in |01). 
These fluctuations lower the energy of the Bogoliubov 
ground state by — iVVb/4 compared to the mean-field 
ground state; the relative reduction in the interaction 
energy is 0(N~ 1 ). The vortex, which becomes ener- 
getically stable [H, [l(| at the critical rotation frequency 
il c = oj± — NVq/AH, is on average slightly off-center by 
0(1/ y/W) (in units of dj_) due to these quantum fluc- 
tuations. In this section, we investigate the nature, in 
real space, of correlations induced by Bogoliubov fluctu- 
ations. 

The Bogoliubov ground state of a single-vortex LLL 
system at ft = fi c is [J] 



|G) 



1 



(5) 



where a m annihilates a particle with angular momentum 
to from the state |0to), and \N\) is a coherent state with 
Ni particles condensed in 1 1 ) , satisfying the eigenvalue 
equation ai\Ni) — y/Ni\Nx). This wave function does 
not conserve the particle number. In order to find its 
form in configuration space, we restrict the number of 
particles to N (assumed to be even) and project |G) 
onto the jV-particle Fock space. This new wave func- 
tion, |G;iV), can be approximated as a sum over states 
with N — 2m particles in |01) and to particles in 1 00) and 
102), 

1 N/2 

|G;A0~-=^ (-l/^) m |TO,7V-2TO,TO), (6) 

v 2 m=0 

where \no,ni,ri2) contains rij particles in |0j) (with j = 
0, 1, 2). The norm of this wave function is (G; N\G; N) = 
1 — 2~' 1+JV / 2 ' and approaches unity when N — > oo. 

The first term in the sum (to = 0) is just the orig- 
inal mean-field many-body ground state, (z|0, N, 0) ~ 
Z\ ■ ■ ■ Zn where z = {zi, Z2, . . . , z^}. The to = 1 term in- 
cludes first-order corrections and yields (z|l, JV — 2, f) ~ 
V \_z\z2 ■ ■ ■ zjv-i(2 2 v r/v / 2)] where V denotes the sum of 
the distinct permutations with respect to the Zj's needed 
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to symmetrize the wave function. After simplifying this 
expression (details in Appendix [B]), we find that the 
first-order Bogoliubov corrections take one pair of parti- 
cles out of the condensate and correlate them through a 
bosonic Jastrow factor, 



where 



Il<l2 k^ix,i2 



(7) 



Similarly, the second-order Bogoliubov correction, the 
rn = 2 term in Eq. ©, leads to (z|2, N -4,2) ~ 
T[z°z%z 3 ■ ■ ■ ZAr_ 2 (zw_i/y2)(zir/V2)] . Simplification of 
the resulting expression (details in Appendix [B]) shows 
that two pairs of particles are correlated through two 
simultaneous Jastrow factors, resulting in the following 
term in the wave function 



z u ) Yl z k = J 2 (z), (8) 



where the primed sum indicates the constraints i\ < 

hi h < U, h < h, «2 7^ U- 

In fact, the real-space projection of the m th term in the 
expansion ^ has up to m simultaneous Jastrow factors. 
To see this structure, we recast this term as 



(z|m, A — 2m, m) 

~ T> [4 ■ ■ ■ z°m(4 l+ i/V2) ■ ■ ■ (zi n /V2)z 2m+1 ■■■zn] 
1 r i ~h 

= ^/2 [(iV-2m)(m l )J OT { 0^,2^,, } ( Z ) > ( 9 ) 



where ( N _ 2m ) ( ) is the number of distinct terms pro- 
duced by the permutations. The monomial symmetric 
polynomial § m a (z) is defined in Appendix |AT and its 
representation in terms of symmetric polynomials with 
successive number of Jastrow factors, determined in Ap- 
pendix [Bj is given by Eq. ([BTO]) . Thus, 



(z|m, A — 2m, m) = 



2 m/2 



AT! 



1 



m! y (A - 2m)! 



(A -2m)! ^ (2m-2i)! 

^ (m-j)\j\ Jj[Z) 



(10) 



where the A- variable symmetric polynomial Jj(z), given 
by Eq. (|A3[) . includes j successive Jastrow factors. We 
immediately sec up to m pairs of Jastrow-correlated 
particles in the m th -order Bogoliubov correction to the 
mean-field ground state. 

Substituting Eq. (| 10[) into Eq. © and changing the 

order of summation using the identity J2m=o £4=0 = 

Ej^o Em=j! we finally arrive at the expansion of the 
Bogoliubov ground state in terms of Jastrow polynomi- 
als, 



(z|G;A) 



N/2 

-^^•(z), (11) 

V2 3=0 



A = 



A! 



2 m m! 

m— 
iV/2 



(N - 2m)! ' 



.4 



(-l) m (2m- 2j)\ /(A- 2m)! 
2 m- J+ i ( m _j)! j! V Wl 



Equation (jllj) shows how incorporating Bogoliubov fluc- 
tuations in the mean-field ground state leads to pairs of 
particles being forced out of the condensate and corre- 
lated in the Jastrow form. The last term in the expan- 
sion above has correlations represented by A/2 Jastrow 
factors, and its coefficient is 0(N~ N ) for large N. 

Note that in the thermodynamic limit (A —¥ oo), 
the mean-field ground state as described by the Gross- 
Pitaevskii equation is the true ground state of the sys- 
tem (see, e.g., Ref. [l|). In fact, the relative reduction in 
the energy between the Bogoliubov and the mean-field 
ground states is 0(N~ 1 ) for large A Q. In mesoscopic 
Bose-condensed systems, the role played by the correla- 
tions can be significant, with the Bogoliubov wave func- 
tion energetically favored over the mean-field solution. 
Moreover, in Eq. (1X1]) . the ratio of coefficients of succes- 
sive terms decreases with increasing j\ most of the re- 
duction in the interaction energy is due to the first term, 
with only a single Jastrow factor [as in Eq. (JT))]. 



IV. EXTENSION TO MEAN-FIELD VORTEX 
LATTICES 

As discussed above, the Bogoliubov ground state |G), 
through the quantum fluctuations, has a lower energy 
than the mean-field ground state. It is clear from the 
form of the Jastrow polynomial Jj(z) in Eq. (|A3|) that 
this lower-energy state |G; N) is constructed by correlat- 
ing j pairs of particles through j distinct Jastrow factors, 
thereby leaving only N—2j particles in the original mean- 
field condensate, |01). As a second example of the effect 
of Jastrow correlations, we argued, using the wave func- 
tions studied in Ref. Q for a gas of attractive bosons, 
that correlating two particles through a Jastrow factor, 
Eq. ([2]), reduces the energy for repulsive bosons. 

We now show that such Jastrow correlations also lower 
the energy of a vortex lattice state. In mean-field theory, 
an A-particle LLL condensate with A v vortices at the 
positions (on a triangular lattice) takes the form 

N N v 

i=l 4=1 

For large N v , the system is well described by the Thomas- 
Fermi approximation (Til [l2j], with the Thomas- Fermi 
radius R and rotation rate fl given by the solution of the 
two equations 



(R/d ± ) 2 = 



I AbNVo 



= (Sl/u ± )N v , (13) 
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where b ~ 1.158 is the Abrikosov lattice parameter. The 
state (T2|) is not an eigenstate of the total angular mo- 
mentum operator L, but has (L) — hN[^(R/d±) 2 — 
1]P. 

To study the effect of Jastrow correlations on the en- 
ergetics of the vortex lattice, we construct a trial wave 
function by removing two particles from the mean-field 
condensate and simutaneously correlating them, arriving 
at the wave function 

ip tl -(z; N v ) = ^2(z il ~z i2 ) 2 tlj nll (z~{z il ,z i2 };N v ) (14) 

il<i2 

where ip m [(z — {z^, Zi 2 }; _/V v ) is an (N — 2)-particle coher- 
ent state (with particles %\ and ii removed) supporting 
the same vortices as the original state (TT2)) . Since the 
Jastrow factors in Eq. (| 14[) force the particles away from 
each other, we expect the cloud for the correlated state 
to extend further in space compared to the mean-field 
one; in fact, the correlated state carrying the same total 
angular momentum as the mean-field one has a radius 
given by 

i? t 2 r ^i? 2 (l + 4/iV). (15) 

The total interaction energy (found after a tedious cal- 
culation, details of which are beyond the scope of this 
paper |13j ) is 

v tr ~V (4bv/3)(N-8), (16) 

where v — N/N v is the filling factor; this result is valid 
for large filling factors. Including Jastrow correlations in 
the trial wave function indeed lowers the energy [albeit 
by a term 0(iV -1 )] compared to mean-field vortex lattice 
state, for which v m { ~ Vo(4bv/3)(N— 1) at the same value 
of the total angular momentum. The relative change in 
the interaction energy is similar to that for a single- vortex 
system as well as that for attractive bosons of Ref. Q. 
Moreover, due to correlations, the average density at the 
vortex cores is non-zero for the trial state (TT4")) , similar 
to the behavior found in Ref. Q for a single vortex. In 
the limit of large number of vortices and for v 3> 1, we 
find that the density at the vortex core is [l3[ ntr(O) ~ 
v~ l |£j| 2 e~'^' (except for the central vortex). 

We note that a relative 0(N~ 1 ) change in the energy 
is not enough, in the thermodynamic limit, to drive the 
system towards the strongly correlated regime where the 
vortex lattice melts Q. A detailed description of the 
melting of the lattice will involve states with large num- 
bers of Jastrow-like correlations, e.g., as in Read-Rezayi 
states. Therefore, vortex lattice wave functions of the 
form pT|) are only good for large filling factors where the 
Gross-Pitaevskii equation is an excellent approximation. 



V. CONCLUSION 



This work is an initial study of the role of correla- 
tions in the ground state of a vortex lattice state, in 
the regime where the Gross-Pitaevskii equation is a good 
first description and quantum fluctuations are small. Al- 
though the advantages of including such interparticle cor- 
relations are clear - keeping the particles apart and re- 
ducing the interaction energy in the system - the de- 
tailed correlations in the exact ground state of the vor- 
tex lattice are not known analytically. Quantum fluctua- 
tions, driving the system towards a melting transition to 
strongly correlated quantum Hall states, become more 
pronounced as the angular momentum per particle ap- 
proaches O(N) and the particle density becomes small, 
underlining the importance of interaction-induced corre- 
lations in this transition. The real-space form of the Bo- 
goliubov ground state of a single- vortex condensate in the 
LLL studied here shows explicitly the Jastrow-like corre- 
lations of pairs of particles in this state. The Bogoliubov 
wave function is a superposition of the original uncorre- 
cted mean-field ground state and correlated states with 
successive number of Jastrow pairs. As we showed, in- 
cluding Jastrow-correlated pairs (similar to those in the 
single-vortex Bogoliubov wave function) in a LLL sys- 
tem with N v vortices on a triangular lattice lowers the 
energy compared to the mean-field wave function with no 
correlations; this state also exhibits non-zero density at 
the vortex cores, reflecting the quantum uncertainty in 
the vortex positions. Generally, interparticle interactions 
lead to the occupation of single-particle states that were 
originally unoccupied in the mean-field picture, allowing 
the system to explore larger regions of phase space, as 
effectively takes place in our trial wave function (fl"4"|) . as 
well as in Ref. [4| in the single-vortex Bogoliubov wave 
function. The next step needed is a systematic study 
of the evolution of the populations of the single-particle 
states of the vortex lattice with increasing angular mo- 
mentum. 
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Appendix A: Monomials and Jastrow factors 



In this Appendix, we define the elementary and monomial symmetric polynomials and find the expansion of the 
latter polynomials in terms of symmetric polynomials with Jastrow factors. We consider a set of N variables, denoted 
by z = {zi, z 2l ■ ■ . , zjv}i and a set of N exponents, denoted by a = a 2 , ■ • ■ , ojat}- The elementary symmetric 
polynomials defined on z are 



5 (z) = 1, Si(z) = ^Zj 1) 5 2 (z) = ^2 Zi t Zi 2 , S 3 (z) = ^2 



Sn{z) 



E 



J1<12<13 



*1<»2<-<«N 



The monomial symmetric polynomials, denoted by tn a {z), are defined as the sum over all z"* 1 z 2 * 



■-N 



k 

where the 



exponents a ix , a i2 , . . . , cti N range over all distinct permutations one can get from a [8| . For example, for N = 3, we 
have m {2 ,o,o}(zi, z 2 , z 3 ) = zfz%z$ + z\z\zl + z?z£zf = z\ + z\ + zf . 

The identity z 2 + z\ = (z\ — Z2) 2 + 2z±z 2 for N — 2 can be rewritten in terms of the symmetric polynomials defined 
above as n^ 2 .Q}(zi, z 2 ) = (zi — z 2 ) 2 + 2 S 2 (zi, z 2 ). There exists a similar identity for N = 4 [see Eq. (|B6[) below] which, 
in the language of symmetric polynomials, becomes 

^{2,2,0,0} (zi,z 2 , z 3 ,z 4 ) =i (zi - z 2 ) 2 (z 3 - z 4 ) 2 + (zi - z 3 ) 2 (z 2 - z 4 ) 2 + (zi - Z4) 2 (z 2 ~ z 3 y 

+ (zi - z 2 ) 2 z 3 z 4 + (zi - z 3 ) 2 z 2 z 4 + (zi - Zi) 2 Z 2 Z 3 

+ {Z 2 - Z 3 ) 2 ZlZi + (z 2 - z 4 ) 2 ziz 3 + (z 3 - z 4 ) 2 z 1 z 2 
+ 6s 4 (zi,z 2 ,z 3 ,z 4 ). (Al) 
We now find a similar identity for general N, assuming, without loss of generality, that N — 2n. Defining 

W{2...2,o...oj(z) = V[z 2 z 2 ■ ■ ■ z 2 n z° n+1 zl +2 ■ ■■z% n _ 1 z% n ], (A2) 



i Jastrow pairs 



Ji(z) = V[(zi - z 2 ) 2 (z 3 - z 4 ) 2 ■ ■ ■ (z 2l _i - z 2l ) 2 z 2i+1 z 2l+2 ■ ■ ■ Z 2n -lZ 2n ] , 



we can write 



m {2 . 



71 

!,0...0}( z ) = Ji(z 

i=Q 



(A3) 



(A4) 



Note that Jo(z) = S 2n (z). 

To find the coefficients, we proceed as follows. First, we set Zj = 1 for all j. Therefore, in the expansion (|A4[) . only 
the Co-term is non-zero. Since the number of terms in the monomial is (2n)!/(n!) 2 and all are equal to 1 in this case, 
we find Co = (2n)!/(n!) 2 . Next, we set z\ = and Zj^i = 1. We find that on the right side of Eq. (|A4[) . only the 
ci-term survives if z\ is one of the two variables in the Jastrow pair, while on the left side, only terms with z\ survive. 
There are ( 2n 1 _1 ) ways on the right to make a Jastrow pair with z\ and one other variable; the remaining variables 
can be arranged in only one way. On the left, for terms with Z®, there are (2n— \)\/\n\(n — 1)!] ways to get a non-zero 
value (which is 1). Therefore c\ = (2n — 2)\/[n\(n — 1)!]. Generalizing this approach to find Cfe (with k < n), we set 
Z\ = z-x = • • • = Zk = and the rest of Zj 's equal to 1 and proceed as before to find 



Ck 



{N-2k)\ 



(2n-2fe)! _ 
nl(n - k)l ~ {N/2)\(N/2-k)\ 



which yields 



N/2 

fft{2...2,0...0}(z) = E 



(N — 2i)\ 



{N/2)\{N/2 i) 



■ Ji(z). 



Let us count the number of terms in each Jj(z). We write the act of the permutation operator V as 

k^i 1 ...i 2 j 



(A5) 



(A6) 



(A7) 
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where the prime on the sum indicates the following conditions 

h < is, h < u, i2j-i<i2j, 

ii < i 3 < i 5 < ■ ■ ■ < faj-i, (A8) 
i2i ^ hl+i, 121+2, ■ ■ ■ > izj for 1 < I < j. 

In order to construct the Jastrow factors, we choose the z's in them as follows. We pick two z's for the first Jastrow 
factor in (^) distinct ways, then the two different z's for the second factor in ( 2~ ) distinct ways, and so on until 
the last one for which there are ( N ~^ +2 ) distinct ways. Therefore, we have N\/[(N — 2j)\2^\ distinct ways to pick 
the z's for the Jastrow factors. Moreover, the Jastrow factors can be permuted in j\ distinct ways among themselves 
while keeping Jj(z) invariant; however, only one of these permutations satisfies the constraints above. The remaining 
z's can be arranged in only one way. Thus, each Jj(z) has 



distinct terms. 



Appendix B: Expansion terms 



In this Appendix, we discuss the method we use to simplify the expansion terms in the Bogoliubov ground state, 
Eq. ([6]), and to bring out the Jastrow factors that include interparticle correlations. The m = 1 term in Eq. (j6]) is 
proportional to V\z\z<i ■ ■ ■ zjv-iz^] where the permutations yield ( N _ 2 ) (J) distinct terms, i.e., 



1/2 



(Bl) 



<*|1,JV-2,1> = [(/ 2 ) ©j V [z° lZ 2 ■■■zn-i (4/^2)] • 
To proceed, we note that the indices of summation (and multiplication) are, in fact, dummy variables and find 

5>m n ^=Ef(*i+^i) n * = e k* - *) 9 + 2 **j] n ( b2 ) 

i^tj k^i,j i^j k^i,j i<j k^i,j 



We thus write 



n 



fj j i -/, 

k 



(B3) 



which leads to the expansion of (z|l, N — 2, 1) in terms of Jastrow polynomials. 

The m = 2 term in Eq. @ is proportional to V\z\z\z-i ■ ■ ■ zn^2Z%-i z n] where the permutations yield ( 2 ) 
distinct terms, i.e., 



<z|2,iV-4,2) 



W-4/ V2 



-1/2 



V [z°z 2 • • ■ z N ^{z 2 N _jV2) (z 2 N /V2)] . 



(B4) 



The permutation operator can be expanded as 



V -z°z°- 



#13^44 



2 2 

2 Z i3 Z ii 



n 



- V ifz°z° 



2 2,0 2 02,0220 

2 Z i3 Z in ' Z ii Z i2 Z i3 Z ii ' Z ii Z i2 Z i3 Z U 



+z 2 z° z° z 2 + z 2 z° z 2 z° + z 2 z 2 z° z° ) 



n 

k^i± . . .24 



z k 



(B5) 



where the unrestricted sum on the left side overcounts each factor of z" ± z" 2 (with a = 0, 2) by 2 (e.g., z^zf; and z^zj); 
as before, the equality originates from the permutations on the dummy variables ii, 13, «4- The terms in parentheses 
above can be rewritten in a more suitable form with the identity 

z\z\ + z\z\ + z\z\ + 2%zl + z%zl + z\z\ = \ \{z\ - z 2 ) 2 (z 3 - z 4 ) 2 + {z\ - z 3 ) 2 (z 2 - z 4 ) 2 + (zi - z 4 ) 2 (z 2 - z 3 ) 2 ] 
+ [(zi - z 2 ) 2 z 3 z 4 + (zi - z 3 ) 2 z 2 z 4 + (zi - z 4 ) 2 z 2 z 3 + (z 2 - z 3 ) 2 ziz 4 + (z 2 - z 4 ) 2 ziz 3 + (z 3 - z 4 ) 2 ziz 2 ] + 6ziz 2 z 3 z 4 . 

(B6) 
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Using this in Eq. (|B5[) leads to (i) 3 equal contributions from the hrst term on the right side of Eq. (IB6|) . each of 
which leads to a factor of 2 x 2 for converting the unrestricted sum to i\ < i 2 and i 3 < i± and another factor of 2 for 
imposing the condition i\ < i^; and (ii) 6 equal contributions from the second term on the right side of Eq. (|B6|) . each 
of which leads to one factor of 2 for converting the sum to i\ < ii and another factor of 2 to count interchangablity 
of Z3 and i^. Therefore, 



£4444 n^ 2 £'K 



k^ii . . .24 

where the prime on the sum indicates the conditions (|A8[) . We now have 



k^ix ...14 *i<*2 kj^ii ,12 k 



z. 



n Zk + e ( z *i ~ z ^) 2 + ( 4) 31 n zk 

k^ii ...14 *i<*2 kj^ii,i2 k 



(B7) 



and, in turn, the expansion of (z|2, N — 4, 2) in terms of Jastrow polynomials. 

As shown in Eq. (0), the m th term in Q is proportional to the monomial OT{ ...o,2...2,i...i}( z ) which we rewrite as 



m. 



•{0...0,2...2,1...1}1 



£ 



fc^l'l ...J2rn 



1 



£ 



(m!) 2 (2m)!/(m!) 2 
(2m- 2j)! 



•{2...2,0...0}( 2: ii • ■ • z i2m) z fc 

fc^il...i2m 



(2m)! z — ^m!(m — j)! 



i 1 ^---^i 2m k^i 1 ...i 2m 



(B8) 



where (i) due to ii, . . . , i2m being dummy variables, we have used a method similar to (|B5|) to get the second equality 
and to represent all the terms in the sum by a new monomial acting on a limited set of z's; and (ii) we use Eq. (|A6[) 
in the last equality. 

To proceed further, we need to recast the square bracket above (which includes Jastrow polynomials defined on the 



subset {zi t , . 



, } C z) in terms of Jastrow polynomials acting on the set z. The result is 



E J i( zi i 



n zk 



(2m) 



i( N 

\2m 



, (2m)!2^ 



(B9) 



where the details of this derivation are as follows. Clearly, the sum over ii, . . . , i2 m leads to an overcounting which we 
need to determine separately for each j. Since i\, . . . ,t2m are dummy variables, each term in Jj(zi 1 . . . Zi 2m ) produces 
the same polynomial after being summed over; this brings in an overcounting factor given by (|A9|) . On the other 
hand, since Y^i 17 n 2 = ^X^^j due to the conditions (|A8|) . we are overcounting by a factor of 2 for each Jastrow 
factor (of which there are j) and by a factor of 2 for each two adjacent Jastrow factors (of which there are j — 1), in 
toto, an overcounting factor of 2 2j_1 . Permutations of the remaining z's outside the Jastrow factors in Jj(zi 1 . . . Zi 2m ) 



leave it invariant, and this leads to an overcounting factor of (2m — 2j)\. Therefore, when we transform Jj(zi 1 



,) 



to Jj(z\ . . . Zff), we overcount by a factor of (2m)! 2 J 1 / j\ for each j ' ^ 0. For the special case of j — 6, since there 
are no Jastrow factors present in Jo(zi x . . . Zi 2 ), we instantly end up with Jq(z\ . . . z^) but overcounted by a factor 
of(2m)!(^). 

Hence, we write the monomial OT{ ...o.2...2,i...i}( z ) m terms of symmetric polynomials with successive number of 
Jastrow factors as 



Tft{0...0,2...2,l...l}( z ) 



— y 



(2m -2j)\ 



m! * — ' (m — 7')! 

3=0 y ■' ' 



N 
2m 



2^ 

Sjo H 77- (1 - S j0 ) 



(BIO) 



and, in turn, find the expansion of (z|m, N — 2m, m) in terms of Jastrow polynomials, Eq. ([T 
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